Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN MfeK style file vl.4) 



Adaptive Optics observations of LBQS 0108+0028: if-band 
detection of the host galaxy of a radio-quiet QSO at z ~ 2 
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ABSTRACT 

We report the first unambiguous detection of the host galaxy of a normal radio-quiet 
QSO at high-rcdshift in X-band. The luminosity of the host comprises about 35% of 
the total i\-band luminosity. Assuming the average colour of QSOs at z ss 2, the host 
would be about 5 to 6 mag brighter than an unevolved galaxy placed at z m 2, 
and 3 to 4 mag brighter than a passively evolved L* galaxy at the same redshift. The 
luminosity of the host galaxy of the QSO would thus overlap with the highest found 
in radio-loud QSOs and radio-galaxies at the same redshift. 
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1. INTRODUCTION 

Recent evidence that the cosmological evolution of the den- 
sity of star formation in the Universe (Madau et al. 1996) 
follows closely the QSO density evolution (Boyle & Terlevich 
1998) emphasizes the need to study the kinds of galaxies that 
host Active Galactic Nuclei in order to understand the link 
between star-formation and nuclear activity, and potentially 
the role of nuclear activity in galaxy formation. 

At the peak value of QSO density (z ~ 2 to 3), the 
few QSO host-galaxies detected so far present rest-frame 
UV fluxes that reach up to 20% of the total QSO luminos- 
ity, indicating star-formation rates about 200 M©/yr and 
above for both radio-loud (Lehnert et al. 1992) and radio- 
quiet samples (Aretxaga, Boyle & Terlevich 1995, Hutchings 
1995). These values are almost an order of magnitude above 
those of field galaxies at similar redshifts selected through 
Lyman Break techniques (Steidel et al. 1996, Lowenthal et 
al. 1997). The properties of these QSO hosts are not unprece- 
dented, since they follow very closely the luminosity-size 
relation of nearby star forming galaxies, overlapping with 
its high- luminosity end (Aretxaga, Terlevich & Boyle 1998) . 
However, the UV fluxes only carry information about the 
high-mass end of the stellar populations in the galaxies, and 
say little about the bulk of the stellar mass which is better 
characterized by optical to NIR observations. 

Although a few hosts of extreme radio-loud QSOs at 
z ~ 2 have been detected in NIR bands (Lehnert et al. 
1992, Carballo et al. 1998), attempts to image the hosts of 
normal radio-quiet QSOs at the same redshifts have been 
unsuccessful to date (Lowenthal et al. 1995, Aretxaga et al. 



1998). Imaging radio-quiet systems, which constitute more 
than 95% of all QSOs, is important in order to character- 
ize the bulk of the population. The observed optical sizes of 
FWHM~ 1 arcsec (Aretxaga et al. 1995), clearly demand a 
technique which offers the highest available angular resolu- 
tion. 

In this paper we focus our attention on the detection 
of the host of a normal radio-quiet z ~ 2 QSO with the 
Adaptive Optics System in operation at the ESO 3.6m tele- 
scope in La Silla. Preliminary results on similar programs 
to image the host-galaxies of QSOs at z ~ 0.5 and 1.7 us- 
ing Adaptive Optics have been presented in a recent confer- 
ence devoted to quasar hosts (Bremer et al. 1997, Hutchings 
1997). 



2. DATA SET: ACQUISITION AND 
REDUCTION 

We selected LBQS 0108+0028 at a(J2000) = l h 10 m 38.1 s 
5(J2000) = 0°44'54" , a V = 18.3 mag QSO at z = 2.005 
which was discovered in the Large Bright Quasar Survey 
(Hewett, Foltz & Chaffee 1995), because it belongs to a nar- 
row redshift-luminosity band (1.8 Ss z <5 2.2, Mb —28 
mag for Ho = 50 Km s _1 Mpc -1 and qo = 0.5) and it lies 
close in projection 9 — 21 arcsec to a bright star of magni- 
tude V — 12.0 mag. The first selection criterion was imposed 
in order to explore the luminosity band that is predicted to 
contain the most luminous hosts by quasar formation theo- 
ries (Terlevich & Boyle 1993, Haenelt & Rees 1993), and it 
has indeed provided a high detection rate of extended fuzz 
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Table 1. Summary of observations 

Object a (J2000) 



LBQS 0108+0028 
star-A 



1 10 38.1 
1 4 1.24 



S (J2000) 

44 54 
55 6.7 



2.005 



V(obj) 
(mag) 
18.3 
13.6 



9(obj - *) 
(arcsec) 
21.0 
20.6 



V(*) 
(mag) 
12.0 
11.3 



t 

00 

10800 
5400 



(Aretxaga et al. 1995). The second condition was imposed in 
order to be able to correct the atmospheric turbulence with 
Adaptive Optics, using a nearby bright reference star on- 
axis, since QSOs at these redshift are too faint to allow for 
direct corrections on themselves. There is no radio-detection 
of this QSO. 

The observations were carried out in if'-band on six 
hours spread through the nights of 1995 October 10th, 11th 
and 12th at the ESO 3.6m telescope in La Silla, with Come- 
On+ (Rigaut F. et al. 1991, Rousset et al. 1994). We used 
Sharp II, the 256x256 Nicmos III array, in the general pur- 
pose 50 marcsec/pixel resolution mode, which gives a 12.8 
x 12.8 arcsec 2 field of view. 

In order to accurately measure a realistic Point Spread 
Function (PSF) we observed every night a comparison dou- 
ble system that mimics the brightness of the reference 
star and its relative distance to the QSO. The compari- 
son system comprises a star of magnitude V = 13.6 at 

a(J2000) = l h 4 m 1.24 s <5(J2000) = 0°55'00.6" , which we 
denote A, separated by 9 — 20.6 arcsec from a star B 
of magnitude V = 11.3 mag at a(J2000) = l h 4 m 1.24 s 

<5(J2000) = 0°55'6.7" . Star B was used as reference star 
to correct for atmospheric turbulence and star A was im- 
aged to serve as a PSF calibrator star in the analysis of the 
QSO profile. 

Coordinates, redshift, l/-band brightnesses, distances 
from reference stars, brightnesses of reference stars, and to- 
tal integration times are summarized in Table 1. 

The observations were carried out in the following se- 
quence: PSF-star (5 x 30 sec) in position 1 - PSF-star (5 
x 30 sec) in position 2 - QSO (10 x 60 sec) in position 1 - 
QSO (10 x 60 sec) in position 2 - PSF-star (5 x 30 sec) in 
position 1 ... in a repeating cycle totaling 3 hr of integration 
for the QSO and 45 min for the PSF-star. Different frames 
were offset by approximately 6 arcsec (distance between po- 
sitions 1 and 2) from one another in order to estimate the 
sky level from contiguous frames. 

The seeing, as recorded by the Differential Image Mo- 
tion Measurement, was 0.8 arcsec during the first night, be- 
ing very stable (±0.03 arcsec), but variable the second and 
third nights (0.8 to 1.8 arcsec). 

The data were reduced with the image processing pack- 
age eclipse^ (Devillard N. 1997). The data was first sky 
subtracted, using an average of contiguous frames with mis- 
placed sources, and then flat-fielded with a gain-corrected 
sky flat frame. Bad pixels were identified in the gain map and 
substituted by linearly interpolated values from nearby pix- 
els. The shift-and-add routines of eclipse were then used to 
register individual frames and coadd them in imaging stacks 



of 10 min for QSO frames and 5 min for PSF calibrator star 
frames. 

The final FWHM of the coadded stacks of the PSF cal- 
ibrator star were 0.3 to 0.4 arcsec during the first night; 0.3 
to 0.9 arcsec during the second night; and 0.7 to 1.0 arcsec 
during the third night. The Strehl ratios attained (8 to 12% 
for Oct. 10th, but below 5% for Oct. 11th and 12th) were only 
acceptable during the first night. The second and third night 
image quality was poor, partly because fast sky variations 
and partly because fast bad seeing, which provoked a reduc- 
tion of the isoplanatic patch to distances much smaller than 
those of our object to reference star systems. 

The images were not flux-calibrated since we didn't ac- 
quire calibration stars and there isn't any A"-band measure- 
ment of this QSO in the literature. 



3. PROFILE ANALYSIS 

For each 10 min QSO images and 5 min PSF-star im- 
ages, we derived an azimuthally averaged radial profile using 
the STSDAS package in irafQ. We checked for variations 
across the chip (position 1 versus position 2) using the PSF 
profiles obtained during the first night, when the seeing was 
stable. We detected no variations within the error bars, and 
therefore coadded contiguous stacks of QSO and PSF-star 
frames, into 20 min and 10 min exposure images. 

A comparison of the QSO profile and the PSF profile, 
normalized in order to reproduce the luminosity of the QSO 
at its centroid, are however different within the inner arc- 
sec. Figure 1 represents two sets of QSO profiles versus the 
PSF-star profiles obtained before and after each QSO ob- 
servation. The recorded profiles were very stable throughout 
the night. 

During the second and third nights the QSO profiles 
are typically indistinguishable from those of the PSF-star 
and, when different, they are enclosed by the varying see- 
ing profiles. Figure 2 shows the comparison between QSO 
and PSF-star for one of these cases. As stated in the pre- 
vious section, during these nights the correction attained 
by Come-On+ was poor due to the poor seeing conditions, 
and there was little improvement in the spatial resolution of 
the images. The resulting stellar FWHM m 1.0 arcsec im- 
plies that we shouldn't expect to have resolved the extended 
structures at r < 1.0 arcsec that we detected during the first 
night, when the seeing was good and stable. 



eclipse is an image processing engine developed at ESO for 
astronomical data reduction purposes in general, and adap- 
tive optics data reduction in p articular (also, see the web site 
tittp : //www. eso . org/eclipse/ ) 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Fundation 
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Figure 1. Radial profiles of two QSO coadded images of 20 min exposure (squares), compared to the stellar profiles (stars) acquired 
before and after the QSO during the first night. Fluxes are normalized to their centroid values 




Figure 2. As Figure 1, for the second night. 



4. RESULTS 

The QSO profile is more extended than that of the PSF-star, 
with K' excesses of about 35% the total luminosity of the 
QSO. This is the first clear detection of the host of a distant 
radio-quiet QSO in the NIR. Previous attempts to detect 
them with direct imaging (Lowenthal et al. 1995, Aretxaga 
et al. 1998) resulted in non-detections at z fa 2, with the 
exception of a marginal detection of a host at r fa 4 arcsec 
(Aretxaga et al. 1998). The upper limits for the luminosities 
set by these studies were an order of magnitude above our 
detection. However, it is clear that this host would have gone 
undetected by these studies, since all the signal is localized 
in the inner 1 arcsec of the QSO profile while, typically, 
the previous studies were carried out with seeing values of 
FWHM between 0.9 and 1.3 arcsec (see also Fig.2). 

If we assume that this QSO has the average colour of 



QSOs at z fa 2 V - K fa 2.2 ± 0.6 mag (Hewett P.C., priv. 
communication) , where the error bar indicates the total am- 
plitude of colours, then the K-h&nd magnitude of the exten- 
sion would be about 17.2±0.6 mag, overlapping with the K- 
band apparent magnitudes of radio-loud QSOs (Lehnert et 
al. 1992, Carballo et al. 1998) and radio-galaxies (Lilly 1989) 
at the same redshift. Our Jf-band detection of the QSO host 
in this analysis demonstrates that there is at least one good 
example of a radio-quiet QSO host with an extremely lumi- 
nous host galaxy at observed NIR wavelengths. 

Since at z fa 2 the observed Tf-band corresponds ap- 
proximately to rest-frame _R-band, we can make an easy 
comparison of our host galaxy with local galaxies observed 
at optical wavelengths. A local L* elliptical has a luminos- 
ity Mft fa —22.8 mag, as derived from the local luminosity 
function of field galaxies (Efstathiou, Ellis & Peterson 1988). 
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Thus, the host of LBQS 0108+0028 is likely to be about 
5 to 6 mag brighter than an unevolved L, galaxy placed 
at z « 2. Taking into account the evolution that the stel- 
lar populations must have experienced between z = 2 and 
z — 0, an L* galaxy at z ~ 2 would be about 2 mag brighter 
in i?-band than nowadays if it had been passively evolv- 
ing since formation (Chariot & Bruzual 1991). The host of 
LBQS 0108+0028 would thus be 3 to 4 mag brighter than a 
passively evolved L* elliptical galaxy placed at z ~ 2. Even 
higher luminosities should be considered if the light we are 
missing near the center of the host is also taken into account. 

The host of LBQS 0108+0028 would also be at least 
5 mag brighter than the average radio-quiet host galax- 
ies of nearby QSOs: < Mv > ~ —21.6 to — 22.6 mag for 
QSOs at < z >Ri 0.2 (e.g. Smith et al. 1986, Hutchings, 
Janson & Neff 1989, Bahcall et al. 1997) which with a 
V — R w 0.7 mag colour for an Sb to E galaxy (Fukigita, 
Shimasaku & Ichiwava 1995) gives rest-frame luminosities of 
M R « -22.2 to - 23.2 mag. Nearby IRAS selected QSOs 
can reside in very luminous galaxies of up to L ~ 6L* (Boyce 
et al. 1996). 

As already noted by Lehnert and co-workers (1992) 
for their radio-loud sample, the density of luminous QSOs 
(Mb Si —28 mag) at z ~ 2 like the one explored in this 
study is about 10 Gpc~ 3 (Boyle et al. 1991), and their hosts 
can be well accommodated at z = by the tail of the lumi- 
nosity function of field galaxies, an idea which has also been 
proposed by Terlevich and Boyle (1993) in their compara- 
tive study of the luminosity functions of QSOs and ellipti- 
cal galaxies. Clearly, a bigger sample of hosts of radio-quiet 
QSOs should be detected in K-band before establishing an 
evolutionary link between these populations of galaxies. 
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